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Introduction 

Since the 1990’s and the COBE experiment, cosmology has entered into an era of high pre¬ 
cision. Measurements of the anisotropies in the Cosmic Microwave Background (CMB) have be¬ 
come increasingly accurate. Observations of Type la supernovae have permitted to put in evidence 
the present acceleration of the Universe expansion. Combined with the distribution of the Large 
Scale Structures (ESS), these observations have lead to the establishment of the standard cosmo¬ 
logical model. Its main ingredients are the dark matter and the dark energy (or a cosmological 
constant), whose origins still have to be understood. But the model faces other theoretical prob¬ 
lems: the CMB temperature is the same in apparently causally disconnected regions in the sky; 
the Universe is very flat and its density must have been extremely fine-tuned to the critical density 
in the past; most Grand Unified Theories (GUT) predict that magnetic monopoles are produced in 
the early Universe and they should dominate the density of the Universe. Moreover the standard 
cosmological model alone does not provide a mechanism for the generation of Gaussian and nearly 
scale-invariant initial density fluctuations. All these problems are naturally solved if the Universe 
underwent an early phase of quasi-exponentially accelerated expansion, called inflation. 

In the inflationary paradigm, Gaussian and nearly scale-invariant density perturbations arise 
naturally from the quantum fluctuations of one (or more than one) scalar field(s) slowly evolving 
along its (their) potential during inflation. Besides theoretical motivations, inflation is supported 
by strong observational evidences. The amplitude Ag and the spectral index Wg of the power spec¬ 
trum of the primordial curvature perturbations have been measured with accuracy by experiments 
probing the Cosmic Microwave Background (CMB) temperature anisotropies and polarization, 
such as the Planck spacecraft Si- the Atacama Cosmology Telescope and the South Pole 
Telescope 0, giving Ag = 2.196('(QQg' x 10^^ and Wg = 0.9603±0.0073 in agreement with many 
inflation models [||]. A strong bound have also been established on the level of local primordial 
non-Gaussianity, = 2.7 ±5.8 S. In 2015, BICEP2 has claimed the detection of primordial B- 
mode polarization of the CMB |^, attributed to the gravitational waves produced during inflation. 
The ratio between tensor and scalar perturbations r = 0.201q qj favors super-Planckian excursions 
of the scalar field responsible for inflation, called the inflaton, and points towards an energy scale 
associated to inflation close to the GUT energy. Since then, it has been argued that galactic dust 
could contribute more importantly to the signal than initially expected, and net detection has been 
transformed in a lower bound on the tensor to scalar ratio, r < 0.11 (95% C.L.). Nevertheless this 
question is still puzzling and future observations are required to affirm or disclaim the discovery of 
primordial gravitational waves iU- 

The simplest way to realize a phase of inflation is to assume that the Universe was filled with 
a scalar field slowly evolving along its potential. Depending on the shape of the potential, the 
primordial scalar and tensor perturbations have different statistical properties. Despite the accurate 
measurements of the scalar power spectrum amplitude and spectral index and the limits on the 
tensor-to-scalar ratio, many single-field models are still consistent with observations, from various 
high energy frameworks. 

In these notes, based on a short series of lectures given in summer 2015 at the 10th Modave 
school of Mathematical physics, I give an overview of the theory of inflation, going from theoret¬ 
ical motivations to the field dynamics, both at the background and linear perturbation levels. A 
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particular care is given to detail the derivation of the scalar and tensor power spectra. Results for a 
simple power-law scalar field potential are presented as a worked example. The second objective 
is to put in evidence and discuss the status of inflation after Planck and BICEP2. After introducing 
basic notions of Bayesian statistics, I will comment the recent results of J. Martin, C. Ringeval, R. 
Trotta and V. Vennin [^, IE] on model comparison. 

The lectures were intended for young researchers working in the domain of mathematical 
physics, not necessary familiar with the theory of inflation. Therefore these notes may be useful 
for master students and young PhD in various fields of cosmology, gravitation, theoretical and 
mathematical physics, wanting to acquire a general and up to date culture on the topic of inflation. 
I tried to make the notes as pedagogical as possible, detailing some derivations, but basic notions 
of cosmology and general relativity are nevertheless obvious pre-requisites. My objective is also to 
give them practical tricks to derive the observational predictions for any single scalar field potential 
that could arise in their work, and to confront their model to the most recent observations. 

After introducing the motivations for an early phase of inflation in Section ||, I will define the 
usual observable quantities in Section ^ In Section ||, the equations governing the homogeneous 
dynamics for a single scalar field model are given and the common slow-roll approximation is 
introduced. The linear theory of cosmological perturbation will be used in Section ^ to derive 
the power spectrum of scalar and tensor perturbations at first order in slow-roll. Observational 
predictions for a simple large field model are calculated in Section ^ as an example. In Section ^ 
a classification of single-field potentials is presented. The most important results of Planck and 
BICEP2 experiments are presented and discussed in Section In Section after an introduction 
to Bayesian analysis for model constraints and comparison, I will comment on the favored single¬ 
field scenarios based on Planck only and on Planck-i-BICEP2. Einally I will mention in Section ^ a 
few open theoretical questions linked to inflation. 


Note that a lot of material has been borrowed directly from the first and second chapters of 
my PhD thesis [[TI|]. These notes are also based on previous but recent lectures notes and reviews 
by D. Eanglois [|I^, A. Einde [13], D. Baumann and E. McAllister [M], to which the interested 


reader may refer to get further information about the topic. An exhaustive analysis of all single 
field inflationary models proposed so far can be found in Ref. As already mentioned, the 
section on model analysis and Bayesian theory is based on Refs. [Ill]- 


Any comments about these lecture notes are welcome. 
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1. Motivations 

The standard hot Big-Bang cosmological model suffers from several problems. In this section, 
they are briefly explained, and then it is shown how a sufficiently long phase of inflation naturally 
solve them. 


1.1 Problems of the standard cosmological model 

1.1.1 Horizon problem 

Before to express and explain the horizon problem, it is useful to define the various notions of 
horizon in cosmology. It is convenient to define the conformal time 



(1.1) 


where a is the scale factor and t denotes the cosmic time. It is the comoving distance covered by 
light between an initial hyper-surface at time t{ and the hyper-surface at time t. Assuming that the 
Universe evolution starts at h, two points separated by a comoving distance larger than the confor¬ 
mal time T] do not have a causal link. Usually, the initial hyper-surface is set at the Planck-time, 
and all points separated by a comoving distance larger than rj are said to be causally disconnected. 
For an observer in 0 at a time io (see Fig. [I]), rj (to) is the comoving radius of the sphere centered in 
O separating particles causally connected to the observer of particles causally disconnected, rj (?) 
is called the comoving horizon or the particle horizon. It is important to distinguish between the 
particle horizon and the event horizon, which is, for the observer, the hypersurface separating the 
universe in two parts, the first one containing events that have been, are or will be observable, the 
second part containing events that will be forever unobservable. Mathematically, the event horizon 
exists only if the integral 



converges. Finally, it is useful to define the comoving Hubble radius, 1 / (aH). It is smaller than the 
conformal time, that is the logarithmic integral of the Hubble radius. 

The horizon problem is linked to the isotropy of the CMB and can be expressed in the following 
way: how to explain that regions in the CMB sky have the same temperature whereas their angular 
separation is too large to correspond to causally connected patches at the time of last scattering, 
assuming that the standard cosmological model is valid down to the Planck time? 

In the standard cosmological model, the early Universe is dominated by the radiation and the 
chemical potentials can be neglected most of the time. One has therefore aT = constant, and in a 
comoving coordinate system, any physical distance grows like 


d{t) = 

The temperature of CMB photons today is To ~ 2.1 K « 2.3 x 10^'^GeV. On the other hand, 
assuming that the expansion rate is governed by the standard cosmological model at every time, 
the radius of the observable universe, i.e. the radius of the spherical volume in principle observable 
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Last scattering surface 



A B Space 


Figure 1: Scheme [ jl^ illustrating the horizon paradox. The CMB is observed from the hypersurface t = tQ. 
The AB' region at last scattering is isothermal in the CMB sky, although it appears to be constituted of 
causally disconnected patches. 


today by an observer at the center of the sphere, is ~ lO^^m. At the time of last scattering 

tLss> the radius of the observable universe was 

dHoithss) ~ 7 X lO^^m . (1.4) 

Under the same assumption, at recombination, the maximal distance between two causally con¬ 
nected points would roughly be 

(^HLssCtLss) ~ 2 X lO^^m . (1.5) 

At last scattering, our observable Universe would therefore have been constituted of about 10^ 
causally disconnected regions. But CMB photons emerging from these regions are observed to 
have all the same temperature, to a 10^^ accuracy. At the Planck time, the number of causally 
disconnected patches would have been much larger, about 10*^^. 


1.1.2 Flatness problem 

The Einstein equations in a homogeneous and isotropic FLRW Universe give the Friedmann- 
Lemaitre (FF) equations 


d 

a 


871 K I ^ 

-^{p + 3P) + U, 
3m^i 3 


(1.6) 

(1.7) 


where // = d/a is the Hubble expansion rate, a the scale factor. A' = 0,±1 is the curvature, p the 
energy density and P the pressure. m^\ is the 4-th dimensional Planck mass, which should not be 
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confused with the reduced Planck mass Mpi = Starting from the Friedmann-Lemaitre 

equations and neglecting the cosmological constant*, one can find the evolution equation for the 
curvature density fix, defined as 

= 0 8 ) 

where Q. = STip/One gefs 

dflx 

^ = (3w + 1)(1-Ok)Ok, (1.9) 

dlna 

where w = P/p is fhe equafion of sfafe paramefer. This equafion is easily infegrafed when w is 
consfanf. One has 

—^^ = (1-Oko)(—) +Oko, (1-10) 

Ok (a) V«o/ 

where Oko is the curvature today. Since it is constrained by observations (|0 — 1| < O.Ol [p^) one 
has roughly at radiation-matter equality 

|O(aeq)-l|<3xl0-^ (1.11) 

and at the Planck time, 

I^^K)-i|<io-'’°. 

If the Universe is not strictly flat, the ACDM model does not 
so small. 

1.1.3 Problem of topological defects 

In Grand Unified Theories (GUT), the standard model of particle physics results from several 
phase transitions induced by the spontaneous breaking of symmetries. Such a symmetry breaking is 
triggered during the early Universe’s evolution due to its expansion and cooling, and they can lead 
to the formation of topological defects like domain walls, cosmic strings and monopoles. These 
defects correspond to configurations localized in space for which the initial symmetry remains 
apparent (see Fig. 

Let us consider the symmetry breaking of a group resulting to an invariance under the sub¬ 
group M': ^ . The vacuum manifold ^ is isomorphic to the quotient group ^/^ 

Domain walls are formed when the Oth-order homotopy group of ^ is not trivial. They can be 
due to the breaking of a Z 2 symmetry, or if the resulting vacuum contains several distinct elements. 
Cosmic strings are formed when the first homotopy group of ^ is not trivial, for instance for the 
breaking scheme U(l) —)• {Id}. Monopoles are formed when the second homotopy group 712 
of the vacuum manifold is not trivial. This is the case for the breaking of a SO{3) symmetry into 
Jif = {Id}. For higher homotopy groups, the resulting topological defects are called textures. 

Groups involved in GUT are such that the first and second homotopy groups are trivial, 
7ri((f) ~ 712((^) Id. In the SM, there remains a U(l) invariance corresponding to electromag¬ 
netism. The first homotopy group of U(l) is TTi \U (1)] ~ Z. Therefore, by using the property of 

'This is a good approximation because A dominates the energy density only at late times. 


( 1 . 12 ) 

explain why the spatial curvature is 
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homotopy groups [ |T^ 


(1.13) 


one obtains that the second homotopy group of the vacuum manifold corresponding to the breaking 
of a GUT group is not trivial. That induces necessarily the formation of monopoles 
However, monopole annihilation has been found to be very slow 


21]. Asa consequence, 


their energy density today should be 15 orders of magnitude larger than the current energy density 
of the universe. Domain walls can also lead to catastrophic scenarios, but they can be avoided in 
the schemes of symmetry breaking in GUT. Cosmic strings are observationally allowed, but their 
contribution to the CMB angular power spectrum [22] is constrained [23]. 




Figure 2: Illustration of the formation of cosmic strings due to the breaking of the group t/( 1) into 
{id}. After the transition, the Higgs field (j) takes a different value at each point in space. When the Higgs 
field makes a complete loop in the field space along a closed path in the real space, there exists a point inside 
the path for which the phase is not defined. At this point, the Higgs field vanished, the symmetry is restored 
and the resulting string configuration contains energy. This process is called the Kibble mechanism (for a 
review, see [p4[]). 


1.1.4 Nearly scale-invariance of the primordial scalar power spectrum 

The density perturbations at the origin of the CMB temperature fluctuations start to oscillate 
when their size becomes smaller than the Hubble radius. On the contrary, the perturbations whose 
wavelength is much larger at recombination have remained constant and thus conserve their initial 
amplitude. In the CMB angular power spectrum, these super-Hubble perturbations correspond to 
temperature fluctuations at low multipoles (I < 20). The CMB temperature fluctuations at large 
angular scales therefore directly probe the initial state of those density perturbations. 
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With CMB observations, it has been established that the primordial power spectrum of density 
perturbations is (nearly) scale invariant. The present measurements of the shape of the primordial 
power spectrum will be given in details in Section This constrains the possible physical processes 
at the origin of the initial density perturbations. 


1.1.5 Absence of iso-curvature modes 


There are two different kinds of primeval fluctuations: the curvature (or adiabatic) and iso¬ 
curvature (or entropic). 

The adiabatic density fluctuations are characterized as fluctuations in the local value of the 
spatial curvature (hence the name of curvature perturbations). By the equivalence principle, all the 
species contribute to the density perturbation and one has for any fluid /, 


5p dtif 5s 

p rif s' 

where s = S/a^ is the entropy density. Furthermore, one can write 

s(^) 0 

\ s J s s^ 


(1.14) 


(1.15) 


That means that the fluctuation in the local number density of any species relative to the entropy 
density vanishes. 

The entropic fluctuations are perturbations for which 5p =0 and therefore they are not char¬ 
acterized by fluctuations in the local curvature (hence the name iso-curvature). They correspond to 
fluctuations in the equation of state. 

With observations one has determined that the CMB temperature fluctuations are sourced by 
curvature perturbations, and one has constrained the possible contribution of iso-curvature pertur¬ 
bations [p^]. The mechanism producing the initial inhomogeneities needs to generate only (or at 
least mostly) curvature perturbations. 


1.1.6 Why are perturbations Gaussian? 

The statistical properties of the CMB anisotropies are encoded in the power spectrum of the 
temperature fluctuations, i.e. in the two-point correlation function in the Fourier space, as well as in 
the three-point (bispectrum), four-point (trispectrum), and higher order correlation functions. But 
if the fluctuations follow a Gaussian statistic, these latter ones are all vanishing. CMB experiments 
have not detected with a high significance a non-zero value neither for the three-point neither for 
higher-order correlation functions. Since the temperature fluctuations in the CMB are induced by 
density perturbations, the mechanism generating them must be to be such that their statistics is 
Gaussian. 


1.2 Cosmic Inflation 

Inflation is a phase of quasi-exponentially accelerated expansion of the Universe. By combin¬ 
ing the Friedmann-Lemaitre equations and assuming K = 0, one obtains a necessary condition for 
inflation to take place. 


d > 0 


P + 3P <0 . 


(1.16) 
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as 


The amount of expansion during inflation is measured in term of the number of e-folds, defined 


A^(t) = ln 


q(0 

Qi 


(1.17) 


where ai is the scale factor at the onset of inflation. 

The inflationary paradigm is well motivated because it provides a solution to all the problems 
of the standard cosmological model that have been mentioned: 


• The horizon problem: Inflation solves naturally this paradox if the number of e-folds of 
expansion is sufficiently large. Indeed, isothermal regions in the CMB apparently causally 
disconnected at recombination, can actually be in causal contact because of a primordial 
phase of inflation. Assuming that the expansion was exponential during inflation. 


a{t) = Uie 


H/St 


(1.18) 


(it will be shown later that this condition is nearly satisfied) one can evaluate the number of 


e-folds required to solve the horizon problem. At the end of inflation, the size of the current 


observable Universe Aho must have been smaller than the size of a causal region at the onset 


of inflation dn,, 


d-y 


(1.19) 


where Uend is the scale factor at the end of inflation. If inflation ends at the Grand Unification 
scale (Pg/d ~ 10^^ GeV), one needs 


A ~ In 


/ TodHoitp) \ 

\TenddHX^i) J 


>57 


( 1 . 20 ) 


where Tp is the photon temperature today, and for which we have assumed dn, (fl) ~ hiTpi/Tend, 
where /pi and Tpi are respectively the Planck length and the Planck temperature. If this condi¬ 
tion is satisfied, the entire observable Universe emerges out of the same causal region before 
the onset of inflation. 


• The flatness problem: During inflation, the Universe can be extremely flattened. Indeed, if 
we assume H to be almost constant during inflation, one has 

|f2K(aend)| = |f2K(ai)|e , (1-21) 

and with A > 70 and a curvature of the order of unity at the Planck scale, the flatness problem 
is naturally solved. 

• Topological defects: During inflation, topological defects are diluted due to the volume 
expansion and are conveniently stretched outside the observable Universe. 

• The primordial power spectrum: Models of inflation generically predict a nearly scale 
invariant power spectrum of curvature perturbations, and thus can provide good initial con¬ 
ditions for the density fluctuations in the radiation era. We will give more details on how to 
calculate the curvature power spectrum in the case of single-field inflation later in the notes. 
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• Gaussian perturbations: In inflation, all the structures in the Universe are seeded by quan¬ 
tum fluctuations. As the Universe grows exponentially, the quantum-size fluctuations be¬ 
come classical, are stretched outside the Hubble radius, and source the CMB temperature 
fluctuations. All the pre-inflationary classical fluctuations are conveniently stretched outside 
the Hubble radius today and can be safely ignored. The Gaussian statistic of the perturbations 
therefore takes its origin in the Gaussian nature of the quantum field fluctuations. 


Iso-curvature modes: Most models of inflation source only curvature perturbations. Never¬ 
theless, for some models (like multi-field models), the iso-curvature mode contribution can 
be potentially important and eventually observable (e.g. in Ref. [^]). In multi-field models, 
these are induced by field fluctuations orthogonal to the field trajectory. 


2. Observables 


The shape of the CMB angular power spectrum is sensitive to the initial conditions of the 
density and curvature fluctuations. Inflation provides these initial conditions and can therefore be 
confronted to CMB observations. 

Observations of CMB temperature anisotropies and polarization have permitted to measure 
precisely the amplitude of the power spectrum of scalar curvature perturbations, its spectral index 
that quantifies the deviation from scale-invariance, and to put strong limits on the ratio between 
curvature and tensor metric perturbations. We skip here the details of the CMB theory relating 
initial conditions to temperature and polarization anisotropies. This is a difficult topic that goes 
beyond the scope of these notes. Instead we focus directly on the quantities that are commonly 
used to constrain inflation. 


2.1 Power spectrum of primordial curvature perturbations 

The primordial power spectrum of the curvature perturbation is defined from, 

(C(k)C(k')) = (27r)3p^(k)53(k-k0 , (2.1) 

where ^(k) is the 3-dimensional Fourier transform of (^(x). The spectral index of this power 
spectrum Us is defined as 


din Ik^Prik)] 
n, = \+ - 


dink 


( 2 . 2 ) 


where k* is a pivot scale in the observable range, k* = 0.05 Mpc^' for Planck. A power spectrum 
increasing on large angular scales (ug < 1) is called red-tilted, if it increases with small scales it is 
called blue-tilted. Deviation from a scale invariant primordial power spectrum have been detected 
by recent CMB experiments. The power spectrum is observed to be red-tilted, and the case «§ = 1 
is today ruled out. The 1-a bound on the spectral index measured by WMAP was [ |l7| ] were 
Ug = 0.968 ± 0.012. Planck has improved by roughly a factor two the measurement of the spectral 
index [||, giving Ug = 0.9603 ± 0.0073. 

On the other hand, the power spectrum amplitude measured by Planck is 


Ag = 


= 2,196 


+0.051 V in^ 
0.06 ^ 


(2.3) 
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2.2 Tensor-to-scalar ratio 

The tensor metric perturbations, characterized by a power spectrum Ph{k) at the end of infla¬ 
tion, can also affect the CMB angular power spectrum, and especially the B-mode polarization of 
the CMB (for details, see e.g. CMB observations by WMAP and Planck have established a 
strong limit on the amplitude of primordial gravitational waves. Usually, this limit is given as an 
upper bound on the ratio r between tensor and scalar power spectrum amplitudes at the pivot scale. 
After WMAP, the 2-a bound on the tensor to scalar ratio was 0 ], 

r= — < 0.24 . (2.4) 

With WMAP polarization data plus Planck observations of temperature anisotropies, this limit is 
reduced to |[l], ^ r < 0.11, whereas BICEP2 claimed in 2014 a tensor to scalar ratio r = 0.201q qj 


2.3 Other observables 

There are other observable quantities of interest to further constrain inflation models. Some of 
these are defined below: 


The running of the spectral index a^ '. it is defined as 

duc 


a, = 


dink 


(2.5) 




Present data are compatible with Wg = 0- Planck la limits give |[T], ^ cUg = —0.0134±0.0090. 

The /nl parameter: this parameter characterizes the amplitude of the so-called local form of 
the bispectrum of (^, 

= ^/nl [Pf (ki)Pf (k 2 ) -b (2 perm.)] , ( 2 . 6 ) 

defined as the Fourier transform of the three-point correlation function. 


nak;) ) = {2nf5^ £k,- • 


(2.7) 


0=1 


o'=l 


A non-zero bispectrum results from non-Gaussian curvature perturbations. Inflation can be 
a source of small non-Gaussianities, but also the reheating phase, eventual cosmic strings, 
and various astrophysical processes. In the squeezed limit, corresponding to ^3 <C ki ~ ^ 2 , it 
has been shown that all single-field models of inflation yield to = ^(1 — Wg ) ~ 0.02 [ | 2 ^ , 


27]. For multi-field models the value can take higher values, potentially detectable by 


experiments. For the other processes, the amplitude should be ffz- ~ 


pIoc 

/NL 


(see [ ]28[ ] for a 

review), thus a convincing detection of ^ 1 would have ruled out most^ single field 
inflation models. For WMAP, the best limit was obtained in Ref. 0 ] 

= 32±21 ( 68 %C.F.) . (2.8) 


The Planck satellite have reduced the error bars by a factor of four and the present limits are 
/‘g= 2.7 ±5.8 ||]. _ 

^However, non-Gaussianities in single field models could be generated by trans-planckian effects [^] or slow-roll 
violation 
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• The Tnl parameter: this parameter characterizes one of the amplitudes of the local-form 
trispectrum of f,, 

Tt; = Tnl [Pf(ki +k 3 )P^(^ 3 )P^(/: 4 ) + (llperm.)] , 
which is the Fourier transform of the four-point correlation function, 

nC(k,)^ = {2nf5^ • 

The present limits are still relatively weak (Tnl < 2800 (95% C.L.)) and in agreement with 

Tnl = 0 


(2.9) 


( 2 . 10 ) 


3. 1-field infiation: Background dynamics 


The easiest realization of the condition (1.16) is to assume that the Universe is filled with an 
unique homogeneous scalar field called fhe inflaton. The Lagrangian reads 






(3.1) 


where V{^) is fhe scalar field pofenfial and g is fhe deferminanf of fhe FLRW mefric. The equafion 
of mofion (e.o.m.) for fhis Lagrangian is fhe Klein-Gordon equafion in an expanding spacetime, 


. dV 

(l) + 3H^ + —=0. 
dtp 

On fhe ofher hand, fhe energy momenfum fensor reads 

2 5^ 


The energy densify and fhe pressure are fherefore 


(3.2) 


(3.3) 


P = y+ 1^(0) , 
P = ^-V{^). 


(3.4) 

(3.5) 


The condition ( 1.16| ) is satisfied if fhe scalar field evolves sufficienlly slowly, so fhaf ^ 
The expansion is governed by fhe Friedmann-Lemaifre equafions 




L 




(3.6) 


d Sk 


a 3m2j 


[-f + V{(P)] 


(3.7) 
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One can rewrite them in e-fold time N 


H^ = 


Sk 


VW 


m: 


VdA^y 


2 ’ 


1 dH 

Hm 


An 

pi 


m 


#v 

dA^y 


1 


d^(j) ^ dtj) "^pidlnV 

yTTd^di^'^ div “ d^ 


o _ I 

Wpi I dN 


(3.8) 

(3.9) 
(3.10) 


and one sees that the scalar field evolves independently of the Hubble rate dynamics. 

3.1 Slow-roll approximation 

For inflation to be very efficient, the kinetic terms in the F.L. equations must be very small 
compared to the potential. The slow-roll approximation consists in neglecting the kinetic terms 
and the second time derivatives of the field, 

In the slow-roll regime, one has therefore 


(^ < 3//0 . 


(3.11) 


H^ = 


871 

3w7p 


VW 


dV 

3Hd = -. 

^ dtj) 


(3.12) 

(3.13) 


(3.14) 


Using the number of e-folds as a time variable, the field evolution is governed by 

d(/) _ 'Wpi 1 dV 
m ~ ""^Fd^ ■ 

One sees that a large number of e-folds is realized in a small range of 0 when the logarithm of the 
potential is very flat. 

The slow-roll regime is an attractor [32] such that typically a few e-folds after the onset of 
inflation, the slow-roll approximation is valid. As shown later, studying inflation in the slow-roll 
regime is very convenient because in this case model observable predictions for the scalar and 
tensor power spectra are directly related to the scalar field potential and its derivatives. Now let us 
introduce the Hubble-flow functions [33], 


H 


d > 0 


t^n+l — 


din |e„| 


dN 


Using these functions, the F.L. and K.G. equations can be rewritten as 


= 


871 V 
m^3 — Si ’ 


^ = 


-1 


dV 


3 + -£2 - £1 


H 


d(/» ’ 


(3.15) 

(3.16) 

(3.17) 

(3.18) 
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and one sees that the slow-roll regime is recovered when 

£i < 3 , £2 < 6 - 2 £i . 


(3.19) 


One sees also that £i < 3 is required for satisfying the condition > 0. In the slow-roll approx¬ 
imation, they can be expressed as a function of the potential and its derivatives. For the first and 
second Hubble-flow functions, one has [ M ] 




m; 


pi 


1 dV\ 


1671 \V d(j) J 




m 




IdVy \d^V 

V^J ~v'^ 


(3.20) 


-h 


The Hubble flow functions are usually referred as the slow-roll parameters. Finally, let remark that 
there exists several definitions for the slow-roll parameters, e.g. £ and T]sr defined as 


£ = £i , 

_ ^ £ "^pi V" 

~ 2He~ 4k V 


(3.21) 

(3.22) 


such that the relation £2 = — fisR + 4£ is verified. 


4. Theory of cosmological perturbations 

A major success of inflation is that it provides a mechanism for the generation of density 
perturbations that seed all the structures in the Universe. The classical density perturbations orig¬ 
inate naturally from quantum fluctuations that grow and become classical due to the exponential 
expansion. 

The theory of cosmological perturbations describes how the scalar field and the metric fluctu¬ 
ations evolve during inflation. At the linear level, the homogeneous metric is perturbed by 5g^y, 

g^jv(x) = + 5g^v{x) . (4.1) 

There are 10 degrees of freedom (d.o.f.) associated to the metric perturbation 5g^v They can 
be decomposed in 4 scalar d.o.f. A,B,C,E, 4 vector d.o.f. B, et resulting from two space-like 
vectors of null divergence, 2 tensor d.o.f. hij resulting from a space-like tensor with vanishing trace 
and divergence. The perturbed metric then can be rewritten as 

ds^ = a^(ri) {-(1 4-2A)dri^ 4-2(diB 4 -Bi)dx‘dri -f [{I +2C)5ij + 2didjE + 2d^iEj) +2hij] dx‘dxj} . 

(4.2) 


4.0.1 The gauge freedom 

One defines a local perturbation of a quantity Q as 

5Q{x,t) = , 


(4.3) 


14 









Inflation after Planck: from theory to observations 


Sebastien Clesse 


where Q{t) denotes this quantity in the un-perturbed homogeneous space-time. Any perturbation 
therefore depends on the choice of the coordinate systems on each manifold. In other words, if 
a coordinate system is fixed for the un-perturbed space-time, one needs to define an isomorphism 
idenfifying the points of same coordinates in the two space-times. The liberty in this choice implies 
that four d.o.f. are non-physical and only linked to the choice of the coordinate systems on the two 
manifolds. 

Let us consider a transformation of the coordinate system 

+ (4.4) 


where is a space-time like vector. can be decomposed in two scalar {T and L) and two vector 


(L,) d.o.f. via 

= = + D'L; = 0, (4.5) 


where D, is defined as the spatial part of the covariant derivative. Fixing this transformation is 
thus equivalent in fixing 4 d.o.f.. Under fhis fransformation of fhe coordinafe system, the metric 
perturbation transforms as 


8giiv 5g^v + ^^gnv, (4.6) 

where is the Lie derivative along ^. The Lie derivative evaluates the change of a tensor field 
along the flow of a given vector field. It is defined as 


p 1 

^ JT^I--- Mp _ lip _ rpli\--- <y... lip ^ pflj _|_ rj^lli... lip 


Vg 


_ ^ rplii... o... lip ^ p^lii j_ 

— S ^cyivi...Vg Z^^vi...Vg ^ ^vi... a... SV;- 


(4.7) 


1=1 


. 7=1 


Applied fo fhe symmefric mefric, the Lie derivative gives 


•^^gjiv — V^ii 1 


(4.8) 


where is the covariant derivative associated with the metric Asa result, one can show that 


scalar, vector and metric perturbations transform like [ jl^ 

A ^ A + r' -b Af'r , (4.9) 

B ^ B-T + L' , (4.10) 

C^C + JifT, (4.11) 

E^E' + L, (4.12) 

E'^E' + V, (4.13) 

B'^B' + V', (4.14) 

hij —)• hij . (4.15) 

In the same way, the perturbation 8Q becomes 

5Q^5Q + ^^Q. (4.16) 


A quantity is called gauge invariant when it is conserved by the transformation of the coordinate 
system, i.e. when its Lie derivative vanishes. Gauge invariants are for instance the Bardeen vari¬ 
ables [p5|] 


d> = A + J^{B-E') + {B-Ey , 
'F = -C-J^{B-E') . 
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If we set r = B — E', L = —E and L- = —B, , one has 


B = E = 0, Bi = 0, (4.19) 

and the scalar metric perturbations are identified with the Bardeen variables 

A = d), 

C = -'B 

This choice is called the longitudinal gauge, which we adopt from now. 

4.0.2 Scalar perturbations 

Once the gauge is fixed, one can study the evolution of cosmological perturbations in the 
linear regime. Interestingly, scalar and metric perturbations decouple and thus one can consider 
them separately. We will show that vector perturbations decay quickly and can be safely neglected. 
Let us consider first the scalar perturbations in the longitudinal gauge. The perturbed metric is of 
the form 

ds^ = a^(ri) [-(1 +2<I>)dT]2 + (1 -2'T)5,-/dbc'dx-^'] . (4.22) 

The scalar field filling the Universe at a given space-time point has an homogeneous part 0 plus a 
small perturbation 50 <C 

0(x,t) = 0(t)+ 50(x,f) . (4.23) 

In the longitudinal gauge, it is identified to the gauge invariant variable 

50g,i, = 5<p + (j)'{B- E') . (4.24) 

After perturbing the energy momentum tensor, the (0,0) and (/,/) first order perturbed Einstein 
equations read 


(4.20) 

(4.21) 


+ + ^ + ^ (4.25) 

"Ip V d0 / 

Aji 

+ = ^0'50 , (4.26) 

'"p 

'T" + 2Jf^'T' + Af'd>'+ d>(2Ar' + Af'2)+ ^v2(<I)-'T) 

= ^ f(t>'5(t)'-(t>'^<P-a^^5<l)^ , (4.27) 

"Ip V d0 y 

where a prime denotes derivative with respect to the conformal time rj , and where = a'la = aEI. 
Moreover, because 57(7 oc 5/ in absence of vector perturbations, one has <I> = Ti. On the other hand, 
the perturbed Klein-Gordon equation reads 

d^V 

50" + 2Ar50'-V250 + a250^ =2(0" + 2Ar0')<I> + 0'(<I>' + 3'T') . (4.28) 
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One sees that 5(p is directly related to <I> and its derivative, so that there remains only one scalar 
d.o.f.. By combining Eq. (|4.25 ) and Eq. ( |4.27 ), by using Eq. ( 4.26 ) as well as the background equa¬ 
tions, and after a bit of algebra, an unique second order evolution equation for scalar perturbations 
can be derived. 


+ 2 





(4.29) 


It is convenient to work in Eourier space, because in the linear regime each mode evolves indepen¬ 
dently and one can solve the evolution equation for each of them. After a Eourier expansion, one 
can define 


jUs = + , 


win 


2 ,2 {as/e) 

(0^ = k^- -—^ , 

ay/e 


(4.30) 

(4.31) 


where k is a comoving Eourier wavenumber, and equation ( 4.291) can be rewritten in a simpler form. 


+ (O^{k,ri)fl, = 0 . 


(4.32) 


This equation is similar to an harmonic oscillator with a varying frequency. Apart in some specific 
cases, it cannot be solved analytically and one has to use numerical techniques. It is also possible 
to solve it analytically after a Taylor expansion at first order in slow-roll parameters. 

Instead of d> or /Ts, it is a common usage to calculate the mode evolution and the power 
spectrum of the curvature perturbation defined as 


Its power spectrum thus reads 




871 ^ 



By using Eq. (4.28D, one can determine that evolves according to 


^ 3(l+w)V^/ ’ 


(4.34) 


(4.35) 


and as long as the modes are super-Hubble fkjM’ <C 1), remains constant in time. Therefore, 
observable modes re-entering into the Hubble radius during the matter dominated era have kept 
the value they had during inflation, when they exit the Hubble radius, independently of the details 
of the reheating phase^ and the transition between inflation and the radiation dominated era. Eor 
1 -field inflationary models, they can be used to probe directly the inflationary era. 

can be identified to the spatial part of the perturbed Ricci scalar in the comoving gauge, in which the fluids have 
a vanishing velocity (ST; = 0). 

^Let notice that a non linear growth of density perturbations during preheating is expected in some models, possibly 
affecting the linear curvature perturbations on very large scales ||^, 
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Inflation Reheating Radiation Matter 



Figure 3: This scheme form Ref. [|3^ illustrates how observable perturbation modes evolve during and 
after inflation. The horizontal axis represents the number of e-folds generated from the onset of inflation. 
Observable modes exit the Hubble radius on a range of about ten e-folds. From this time, inflation still lasts 
from 50 to 70 e-folds, depending on the energy scale of inflation and on the duration of the reheating 
phase. The curvature and tensor perturbations are constant for super-Hubble wavelengths, until they re-enter 
into the Hubble radius during the matter/radiation dominated era. 


4.0.3 Tensor perturbations 

The metric for the tensor perturbations reads 

As^ = a^{r\) [—dt]^ + (1 +hij)Ax'Ax^^^ , (4.36) 

and the metric perturbation hij is gauge invariant. It is convenient to express the two d.o.f. in /i,y as 

hx o\ 


^ij — ^ 


hx h+ 0 
0 0 0 


(4.37) 


As for the scalar perturbations, one can then write the first order perturbed Einstein equations, 


h''+2J^^ha+V^ha=0, 


(4.38) 


where a = +, x. By defining 

after Fourier expansion, these two equations reduce to 


fl'' + (Otik,T])iJ.i = 0 , 


(4.39) 


(4.40) 
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where 


(of{k,ri) = k^ - 


(4.41) 


The variable h = hijd'f is the analogous of ^ for the tensor perturbations and has similar properties. 
Its power spectrum reads 


2k^ 

<^u{k) = ^ 




(4.42) 


4.0.4 Vector perturbations 

The metric for the vector perturbations in the longitudinal gauge reads 

ds^ = a^{7i)[—dTfi^+ 2d(^iEj-^dx‘dx-’] , (4.43) 

and the vector perturbations Ej can be identified in this gauge to the gauge invariant variable 

^i=Ei-Bi. (4.44) 


The perturbed energy-momentum tensor for a scalar field does not contain any source of vector 
perturbations and the first-order perturbed Einstein equations read 

= 0 . (4.45) 

Vector perturbations therefore decay quickly, since <!>• oc and because a grows nearly exponen¬ 
tially with the cosmic time. That is why vector perturbations are usually neglected. 


4.0.5 Quantification of perturbations 

In the context of inflation, quantum fluctuations are responsible for large scale structures of 
the Universe observed today. The canonical commutation relations are the basis of the quantization 
process. But to define them, one needs the canonical momenta, and thus the action. It is incorrect to 
interpret directly the classical equations of motion [Eqs. ( |4.32 ) and ( 1.40 )] quantum mechanically. 


because it leads in general to an incorrect normalization of the modes [40]. 


Scalar perturbations: If we perturb the total action of the system up to the second order in the 


metric and scalar field perturbations, one finds [ 10 j 




d!^x 


./\2 


2 


-d'^diVdjV + -^V 

Zs 


(4.46) 


where v = a{5(pg± + can be identified to in the longitudinal gauge, and is 

the so-called Mukhanov-Sasaki variable. The quantity Zs is defined as Zg = V^a(j)'jJlf = a^^. 
As expected, the e.o.m for this Eagrangian reads 

v" - -b V = 0 . (4.47) 


The first step of the quantization process is to determine Ti, the conjugate of v. 


% = 


5v' 


= V 


(4.48) 
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Then the Hamiltonian reads 

H = Jdx‘^(^n^ + d‘JdivdjV - . (4.49) 

In a quantum description, the classical variables v and 7i are promoted as quantum operators v and 
71, satisfying the commutation relations 


[v(x,T]),v(y,T])] = [7f(x,T]),7f(y,T])] = 0 , (4.50) 

[v(x,T]),jf(y,T])] =/5(^)(x-y) . (4.51) 

In the Heisenberg picture, the operator v can be expanded over a complete orthonormal basis of the 
solution of the field equation Eq. (1.47). If one takes a basis of plane waves, one has 

1 


v(x,T 7 ) = 

and the equation for the Vi(T]) is 


( 271 ) 3/2 


d^k (v^e^'^'^dk + v|e , 


(4.52) 


4'(T7)+(^"-y)v, = 0. (4.53) 

If the normalization condition 

vU^K(r])-vr(T])v,(T]) = 2/ (4.54) 

is satisfied, fhe creafion and annihilafion operators dk and d^ satisfy fhe sfandard commutation 
relations 

[dk,dk>] = [d+,d+] = 0 , [dk,d+] = 5(3)(k-k') . (4.55) 

At a time rji, the vacuum |0) can now be defined, such fhaf for all k one has 


dk|0)=0. 


From Eq. ( |4.53 ), in fhe sub-Hubble regime, we have 


lim Vk{ri) 

k/aH^-\-oo 


Q-ik{ri-r],) 

Vlk 


This can be used to give consisfenf initial conditions to Eq. (4.32|) 


(4.56) 


(4.57) 


Tensor perturbations: The quantification of the tensor perturbations is analogous. One can first 
determine the second order perturbed action 







a=+,x • 



S'-’dihadjha + —h 
a 


2 

a 


(4.58) 


The perturbations /ice(t],x) are the canonical variables. They are promoted as quantum operators 
and are expanded in plane waves, 

T]) = ^2n)^/2 j. (4.59) 
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The e.o.m. are 


^k.y + ( -] — 0 , 


(4.60) 


similar to Eq. ( |4.40| ). The quantification process can be used to determine the sub-Hubble tensor 
perturbation evolution, 

Q-ikfl-nl) 


lim hkj{ri) = 


y/2k 


(4.61) 


4.0.6 Expansion in slow-roll parameters 


Eqs. (4.32) and (4.4C) can be solved analytically if these are expanded at first order in the 
Hubble flow-functions around some pivot scale. To do so, let us first rewrite 


V 



da _ 1 r £i 

Ha2~ aH^J a^H' 


(4.62) 


In the slow-roll approximation, one has |£, | <C 1. By definition, the derivative of the first and second 
Hubble-flow functions with respect to the number of e-folds are second order in |£,|, 


d£i 

div 


— £i£2 , 


d£2 

m 


= £2^3 • 


(4.63) 


One can therefore neglect their variation over the time taken for observable modes to exit the 
Hubble radius (it corresponds typically to AN ~ 10 [|^). In this approximation, and by using 
Eq. (Q, one thus has 


rr 1 klH 

aH —-1- 

V V 


£i 

a^H 


, 1 aH 

da = -1- 

V V 


£l , 1 , rr —1 -|- £l 

—dt ~-h aHe\ ~-. 

a 7] T] 


(4.64) 


By integrating the last expression, the scale factor is found to behave like 


a(T])~/o|T]| 




(4.65) 


where Zq is an arbitrary parameter. Instead of choosing an arbitrary scale, it is more convenient to 
chose an arbitrary conformal time T]*, and to relate it to the scale Zq via 




1 -|- £i* 

at]* 


^ r[l + ei*(l+ln|T]*|)] , 

Zo 


(4.66) 


where a star subscript denotes the evaluation at the time q*. We have fixed q* such that the follow¬ 
ing relation is verified 

k* = a(T]*)//(T]*) , (4.67) 

where is the comoving pivot mode introduced in Section ^T| . 


The scalar and tensor perturbations evolve according to Eq. ( |4.32 ) and Eq. ( 4.40| ). These 
equations can now be expanded at first order in slow-roll parameters and then solved analytically. 
The first step is to use Eq. (|4.64 ) to expand 


{as/s)" 2 -|- 3 £i -|- ^£2 

as/s ~ t ]2 

a^^ 2 -\- 3£i 

_ _ 

a ~ 


(4.68) 

(4.69) 
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In the approximation that the slow-roll parameters are constant in time, a general solution to ( |4.32 ) 
and ( 4.4C ) can be found 


jUs,t(^T]) = > (4-70) 

where Vg = — | — £i — ^£2 and Vt = — | — £ 1 . It is convenient to express the Bessel function Jv{kr]) 
in terms of the Hanckel functions of first and second kind {kr\ ) et {kr\ ). The quantification 

of the perturbations provide the initial conditions. By using the asymptotic behavior of the Hanckel 
functions, 


V nz 

Hi^\z^oo) = ^ ( 4 _ 72 ) 

V Ttz 

(4.73) 


and by comparing with the Eqs. (4.57) and (4.61), A and B can be determined. For scalar perturba¬ 
tions, one has 


A = 2/—-\/ksin(7rVi)e'(2''' 4+^7.) ^ 

nipi 

B = . 


On the other hand, one can use the limit condition 


r.{l) 

7^1/2-v 

as well as the recurrence relation 


(,^o) = --r --V - 


n 


Z\ v-4 

'2 


r(z+£) = £t/r(z)r(z)+r(z), 

Vr(l/2) = -yEuier-21n2 , 


(4.74) 

(4.75) 

(4.76) 


(4.77) 

(4.78) 


with Thuier — 0.5772 and where t/r(z) is the polygamma function, to obtain the super-Hubble be¬ 
havior of the perturbation modes. Since observable modes are super-Hubble at the end of inflation, 
one obtain the power spectrum expanded at first order in slow-roll parameters around q*, by using 
Eqs. ( [4.65 ) and (4.66). For scalar perturbations, one obtains 

2 


^f(k) = 


Ms 


Hi 


nm^Ei 


1 — 2(C + 2)£i^, + C£ 2 * — (2£i* -|- £ 2 *) lu 


For tensor perturbations. 


2k^ 

^,{k) = ^ 


n 

167/2 


Km: 


P L 


1-2(C+l)£i*-2£i*ln ( — 


(4.79) 

(4.80) 

(4.81) 

(4.82) 
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with C = "Shuler + 2 In 2 — 2. All the slow-roll parameters are evaluated at T]*. At first order in 
slow-roll parameters, the scalar spectral index is therefore 


Ug 1 - 2£|;): £ 2 * ■ 


(4.83) 


For the tensor perturbations, it is 

fit = . (4.84) 

Finally the ratio between the tensor and scalar power spectrum is given by 

r=l6eu. (4.85) 

The amplitude and the spectral tilt of the scalar and tensor power spectra can thus be derived 
easily in the slow-roll approximation, at first order in slow-roll parameters, for a given scalar field 
potential. They are given only in terms of the potential and its derivatives with respect to the scalar 
field. 

Finally, note that r = —Sui is a generic prediction of single field inflation. Checking this 
relation will be a major goal for future experiments. If observations show that it is satisfied, this 
should be seen as a proof that inflation really took place, since most of the alternatives predict 
different behaviors. This relation is called the consistency relation of inflation. 


5. Worked example: the power-law large field potential 

To illustrate the results of the two previous sections, let consider one of the simplest potential, 
of the power-law form (the model is often referred as large field inflation or chaotic inflation): 




(5.1) 


The background dynamics in the slow-roll approximation is given by Eqs. ( |3.12 ) and ( |3.13 ). The 
number of e-folds realized from an initial field value (/>i can be determined analytically. 


N{^) = 


2p 


A 

M„ 


A 

M, 


(5.2) 


The first and second slow-roll parameters read 


p^Ml 


(5.3) 


2pMl 


(5.4) 


Inflation stops when the first slow-roll parameter reaches Ci = 1. This corresponds to the inflaton 
value 

0end _ P 

^“7! ■ 


(5.5) 
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For a given number of e-folds N^, between the Hubble exit of the pivot mode and the end of inflation, 
the inflaton value 0* and the slow-roll parameters can be obtained in a straightforward way. For 
= 60 and p = 2, they read 

= \I^P + , (5-6) 

eu — 0.0083 , £2* —0.0166 . (5.7) 

It is then straightforward to derive the scalar power spectrum spectral index and the scalar to tensor 
ratio, 

Us = 1 - 2 £i* - £2* ^ 0.967 , r = 16 £i* ~ 0.13 . (5.8) 

These predictions are independent of the mass of the field and correspond to a point in the (ug > r) 
plane. Nevertheless, they depend on the reheating history through N*. The mass scale is fixed by 
the scalar power spectrum amplitude given in section One gets M ~ lO^^nipi. For large field 
models, inflation takes therefore place close to the GUT scale. Let remark that the inflaton field 
must be initially super-Planckian in order for inflation to last at least 60 e-folds. Nevertheless, the 
energy density remains much smaller than the Planck scale. General Relativity is thus valid and no 
effect of quantum gravity is expected. 





0/'Mp 



Figure 4: Top: Scalar field potential and its logarithm for the large field model and for p = 2. Bottom-left: 
evolution of slow-roll parameters £i (solid line) and £2 (dashed line). Inflation stops when £1 = 1 (dotted 
line). Bottom-right: evolution of 0(N), for an initial field value = 5mpi 
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6. Classification of single-field potentials 


There exists a large variety of single-field models of inflation, from the simplest power-law 
potentials to more complicated potentials arising in various high energy frameworks. For an ex¬ 
haustive list of models and their analysis, one can refer to the Encyclopaedia Inflationaris [0]. 
Note that all those models are implemented numerically within the ASPIC library. 

In this section, I introduce a classification of single-field potentials, depending on the values 
taken by the scalar field and the slow-roll parameters. It is closed to the Schwartz-Terrero-Escalante 
classification 011111. Note however that it is not really complete since more complicated forms 
of the potential could be envisaged. Moreover, some scalar field potentials can belong to several 
classes depending on their parameters. 

1. Small-field models: Slow-roll parameters are usually such that 0 < £i ^ £ 2 ^ meaning that 
the kinetic energy increases, as well as the ratio between the kinetic energy and the total 
energy. The Higgs model as well as plateau potentials, e.g. of the form F oc [1 — (0/Mpi)^], 
belong to this class. The low value of £i ensures that many e-folds can be realized close 
to ^ = 0 at sub-Planckian field values. The spectral index deviates from unity due to the 
concavity of the potential, which is negative (since £2 ~ —V”jV). 

2. Large-field models: In this class of models, one has usually 2£i > £2 > 0, meaning that 
the ratio between kinetic energy and the total energy increases whereas the kinetic energy 
decreases. The power-law potential (V (j)^) belongs to this class. The spectral index devi¬ 
ates from unity mostly due to the £1 parameter, and given the present constraints this implies 
super-Planckian field values. Note however that the energy density during inflation is still 
much lower than the Planck scale. 


3. Hybrid models: In this class of models, inflation takes place in a false vacuum dominated 
regime, often along a nearly flat valley of some multi-field potential. Inflation can end due to 
the presence of auxiliary fields inducing a tachyonic waterfall instability below some critical 
field value (pc, or due to a change in the shape of the potential. Both the kinetic energy and 
the ratio between kinetic and total energy decrease. The potential is e.g. of the form V 
[I + {^/IJ.)P] (original hybrid model) and many e-folds of inflation can be realized in the false 
vacuum, at sub-Planckian field values. However in this model one gets £2 < 0 and £1 <C 1, 
and thus the spectral index is blue, which is now ruled out. There are nevertheless hybrid 
models, e.g. SUSY models (F-/D-term models and variants), where a flat direction of the 
potential is lifted up by logarithmic radiative corrections. This leads to £2 > 0 and thus a red 
spectrum possibly in agreement with observations. Finally note that in hybrid models, multi- 
field effects can be important, e.g. before field trajectories reach the inflationary valley [^, 
141, or during the final waterfall phase |^, Ih, £7]. 


7. Observational constraints from Planck and BICEP2 
7.1 Planck results 

In 2013 the Planck mission has delivered its measurements of the CMB temperature anisotropies. 
The CMB angular power spectrum has been measured with a high accuracy down to very small 


25 












Inflation after Planck: from theory to observations 


Sebastien Clesse 


scales. It is found to be in a very good agreement with the best fit of the standard cosmological 
model, assuming that primordial perturbations are well described by the amplitude and the spec¬ 
tral index of their power spectrum. With Planck combined to observations of the Baryon Acoustic 
Oscillations and type 1-a supernovae, the cosmological parameters have been measured with an 
unprecedented accuracy ||T], As already mentioned, the scalar power spectrum amplitude and 
spectral index are respectively given by Ag = 2.196^QQg' x 10^^ and = 0.9603 ±0.0073. An 
important point is that < 1 at more than 4a, which rules out models predicting a scale invariant 
or a blue tilted scalar power spectrum (such as the original hybrid model). An upper bound on 
the tensor to scalar ratio (r < 0.11) has also been derived by combining Planck data to the mea¬ 
surements by WMAP of the CMB polarization. The Planck limits in the plane {n^,r) are given in 
Fig. 1^ together with the predictions of some of the most well-known inflation models. One can 
observe that the limit where convex potentials will be disfavored at the 95% C.L. is not far, and that 
simple potentials like V oc (j)^ and V oc (j)^ are already strongly disfavored. Simple supersymmetric 
models like F-term and D-term inflation are also in strong tension with Planck data, because they 
predict typically that 0.98 < Us ^ 1- 

Planck also give the best constraints on the level of local non-Gaussianities, with a bound 

= 2.7 ±5.8 1^. This excludes many multi-field inflation scenarios and tends to favor single¬ 
field models. Planck has also improved the limits on a possible contribution of a cosmic string 
network to the CMB temperature anisotropies. 

We review in more details the implications of the Planck results for inflation models in the 
next section. 


Angular scale 


90° 18° 1° 0.2° 0.1° 0.07° 



Figure 5: CMB temperature anisotropy angular power spectrum seen by Planck [|I|, with the predictions 
for the best fit of the standard cosmological model parameters, assuming a scalar power spectrum described 
by an amplitude As and a spectral index «s. 


26 









Inflation after Planck: from theory to observations 


Sebastien Clesse 



0.94 




0.96 

Primordial tilt (f7s) 


0.98 


1.00 


P/anc/c+WP+BAO 
P/anc/c+WP+h igh L 
P/anc/t+WP 
Natural Inflation 
Hilltop quartic model 
Power law inflation 
Low scale SSB SUSY 
Inflation 

V <x 

V oc 

V !X<j> 

V IX 
A/.=50 
A/, =60 


Figure 6: l-ci and 2-a contours in the plane {ns,r) from Planck plus BAO or WMAP polarization [|l 
The predictions for several inflation models are also reported. 


7.2 BICEP2 results 

The BICEP2 experiment is based at the South Pole in Antarctica and measures the CMB polar¬ 
ization on large scales, in a small patch of the sky. In March 2014 it has reported the measurement 
of B-mode polarization of the Cosmic Microwave Background 0 and claimed that it is of pri¬ 
mordial origin. The signal can be attributed to the tensor modes generated during inflation for a 
tensor to scalar ratio r = 0.201q qj. They also claimed to exclude r = 0 at more than 7a. Combined 
Planck and BICEP2 constraints in the plane are reported on Eig. ^ Both experiments are in 
agreement but one can note some tension between them. As explained in the next section, this is 
translated in many models of inflation that are incompatible with both Planck and BICEP2. 

The theory of the CMB polarization is rather complex and goes beyond the scope of these 
lecture notes. Here I only briefly comment the principal result. The polarization of the CMB is 
due to the Compton diffusion that tends to polarize the radiation in the orthogonal direction to the 
diffusion plane. In the case of perfect isotropy of the radiation, or in the case of dipolar anisotropy, 
there is no net effect and the CMB would not be polarized. However, for quadripolar anisotropies, 
the CMB radiation becomes polarized in averaged. On can distinguish two types of polarization, 
designated by E-mode and B-mode by analogy to the electromagnetism because they correspond 
to curl free gradient and divergence free curl polarizations. The important point is that tensor 
perturbations generate both E-mode and B-mode polarization whereas scalar perturbations generate 
only E-modes. By measuring primordial B-modes, one therefore has access to the amplitude of the 
power spectrum of tensor perturbations, which has been calculated in a previous section. 

Measuring the tensor to scalar ratio has profound implications for inflation. Indeed, the tensor 
to scalar ratio is directly related to the first slow-roll parameter, r = 16ei. Knowing the value of 
Cl* together with the amplitude of the power spectrum of scalar perturbations then permits to fix 
as well as the value of the field and of the potential at the time of Hubble exit of the pivot 
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Figure 7: Planck and BICEP2 combined constraints in the plane {ns,r). Figure from [^. 


scale. Measuring r therefore gives direct access to the energy scale of inflation. Combined with 
measurements of the scalar spectral index, this gives also the curvature of the potential via the 
second slow-roll parameter £ 2 - The energy scale of inflation is given by 



(7.1) 


and therefore BICEP2 results point toward an energy scale of inflation close to the GUT scale. One 
can also note that BICEP2 results lead to a relatively large value of £1 which implies superplanckian 
field excursions and excludes the class of small field models. 

7.3 Controversy about BICEP2 results 

It must be noted that BICEP2 results are still controversial. Planck dust polarization maps 
have been recently published, and they suggest that the polarization fraction has considerable un¬ 
certainties, higher than the previous predictions that where used by BICEP2. In the case there is no 
assumption about dust polarization, except the power spectrum shape 


/^BB ,(iiist 1 —2 3 

C, cx: / 


(7.2) 


then it has been shown that solutions with no gravitational waves and r < 0.11 are favored at 2a 
level. It is currently impossible to affirm with certainty that gravitational waves from inflation have 
been measured One has to wait early 2015 for the Planck results on the CMB polarization. 

Planck is able to detect a tensor to scalar ratio r ~ 0.2, but not r ~ 0.1. Therefore it is likely that 
we will have to wait for further observations, e.g. by BICEP3. Note also that a joint more precise 
analysis between Planck and BICEP has been planned. 

7.4 Future experimeuts 

There exists several projects of CMB experiments that could take place in the next 15 years. 
Besides BICEP3, one can mention EiteBIRD [^], the Cosmic Origin Explorer (COrE) p9[ ] and the 
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Polarized Radiation Imaging and Spectroscopy (PRISM) missions [50], the latter two being now 
gathered in the COrE+ proposal. One aim is to detect primordial gravitational waves at the r ~ 10^^ 
level. It is also possible to have access to a broader range of wavelength modes by measuring 
spectral distortions of the CMB black-body spectrum, e.g. with the Primordial Inflation Explorer 
(PIXIE) [^] or with PRISM. Silk damped acoustic oscillations induce some energy injection in 
the CMB monopole, which results in spectral distortions that can therefore probe the scalar power 
spectrum on much smaller scales than with CMB temperature anisotropies. A variety of inflation 
models could be tested with CMB distortions [ |5^ 

Prom 2020, additional constraints should come from large scale structure observations with 
Euclid, which is expected to improve by a factor 2-3 the present limits on the parameters describ¬ 
ing the shape of the scalar power spectrum [53]. The 21cm signal from reionization and from the 
end of the dark ages could be also a powerful probe for cosmology (see e.g. Refs. [M, 56]), 

and especially for inflation. Puture giant radio-telescopes like the Square Kilometer Array (SKA) 
should inaugurate the detection of the 21cm signal, and an interesting concept of full digital kilo¬ 
meter size radio telescope dedicated to 21cm cosmology has been proposed [^], the Past Pourier 
Transform Telescope (PPTT). 


8. Model comparison: a Bayesian approach 

Considering only single field models, a plethora of scalar field potentials have been proposed, 
arising in various high energy frameworks like axions, non-minimally coupled Higgs field, super- 
symmetry, supergravity. Grand Unified Theories, extra-dimensions, string theory, brane cosmol¬ 
ogy, loop quantum gravity. Most models proposed so far have been listed in Ref. [0 and their 
compatibility with observations has been analyzed. But in order to compare models and to hunt 
what is the best one given the data, one needs to use Bayesian statistical methods. 

In this section, we briefly review the basics of Bayesian inference and define the Bayes factor 
that gives the posterior odds of some model given a reference model. Then, the results for single 
field models of inflation and Planck data are presented. Even if they must be taken with caution, 
the results for a joint Planck and BICEP2 analysis are also discussed. 

This section is mostly based on the recent work by J. Martin, C. Ringeval, R. Trotta and V. 
Vennin [||, |T^, to which we invite the interested reader to refer for further details. 

8.1 Notions of Bayesian inference 

Bayesian inference is based on the Bayes’ theorem that determines how likely are some hy¬ 
pothesis given new data and some prior information. To illustrate how works the theorem, let first 
consider a simple wikipedia example: suppose that a friend told you he met someone at a party. 
If you don’t have any further information and assuming that your friend does not talk more likely 
to women than men, then the probability that the person is a women is 50%, which is denoted 
P(VT) = 0.5 and P{M) = 0.5 (M for man and W for woman). Now if you understand in the con¬ 
versation that this person had long hair, you may want to calculate the probability that the person 
was a woman. If it is known that 75% of women have long hair, which we denote P{L\W) = 0.75 
(E for long hair), and that only 15% of men have long hair, P(L|M) = 0.15 (those are conditional 
probabilities), then you can calculate the probability that a person is a women given that he or she 
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is a long-haired person, P{W\L). It is obvious that P{W\L)P{L) = P{L\W)P{W), which can be 
rewritten 


PiW\L) = = 0-75 X 0.50 ^ 5 

^ ' ’ P{L\W)+P{L\M) 0.75x0.50 + 0.15 x0.50 6 


(8.1) 


Now let us consider a model M, having continuous parameters Oij, and that we want to evaluate the 
probability of some parameter values given some new data D. The Bayes’ theorem is generalized 
as follows: 






( 8 . 2 ) 


where = P{D\dij,^i) is the the so-called likelihood function for the model parameters. 

The function K contains the prior information on the model parameters. S is called the Bayesian 
evidence and is defined as 


^(D|M,-) = J deij^{dij)7iieij\^i). (8.3) 

Note that it is just a normalization factor, and therefore it is not required to evaluate the Bayesian 
evidence if one just want to constrain the parameters dij. On the other hand, if one needs the 
posterior probability of some model, one has to evaluate 






(8.4) 


and then the Bayesian evidence needs to be computed. In one consider that all single field models 
are known and fhaf none of fhem is a priori favored, one has = 1 /Nmodeis- The posterior odds 

of some model ^ compared fo anofher reference model ^Ref are encoded in fhe Bayes factor 

p{jr^,i\Dy ^ ^ 

If is inferesfing fo nofe fhaf fhe Bayes facfor fake info accounf fhe Occam’s razor effecf, in fhe 
sense fhaf models wifh less parameters will be favored againsf models wifh more paramefers and 
same predictions. Anofher imporfanf poinf is fhaf fhe Bayes facfor depends on fhe prior given fo 
fhe paramefers Bij. Differenf priors change fhe model Bayesian evidence. Thus for inllalion, if 
is acfually possible fhaf a same scalar field potential arising in differenf frameworks, mofivafing 
differenf priors, has differenf Bayes facfor depending on fhe considered scenario. 

The principal difficully is fo estimate fhe likelihood function ^{Bij) = P{D\dij,^i) in fhe 
full parameter space, for each model. Usually, fhis is done by using Markov-Chain-Monfe-Carlo 
(MCMC) mefhods. Compared fo sfandard Monle-Carlo fhaf samples fhe par'amefer space wifh 
poinfs randomly disfribufed, a Markov chain in which each poinf depends on fhe previous one is 
builf and fhe sfafisfical disfribufion of fhe poinfs converges fhrough fhe likelihood function. The 
main advanfage of MCMC mefhods is fhaf fhe convergence is nearly linear, which allows fo probe 
high-dimensional parameter spaces. 


8.2 The best inflationary model after Planck 

The Bayes actor has been calculated in for Planck data and for all the 193 single field 
scenarios listed in Ref. 10], using the Higgs model (the only single field model with no free 
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parameter) as a reference. Those results are reported in Fig. One can use the Jeffrey’s scale 
for evaluating the Bayesian evidence between two models: if | lnB[^gj-| < 1 the result is said to be 
inconclusive^ m the range 1 ^ <2.5 there is a weak evidence', for 2.5 < | InB^gjI < 5 there 

is a moderate evidence', if | InSJ^g^l > 5 there is a strong evidence. 

The first important result is that Planck data rule out at a strong evidence level about one third 
of the models. This demonstrates that Planck results are very impressive. Only one fourth of models 
are found to be inconclusive compared to the best model. Interestingly, all the favored scenarios 
have a scalar field potenfial of plateau-type. Among them, Higgs inflation (or equivalently the 
Starobinsky model) is at the top, whereas it is totally predictive since it has no free parameter. 

In Higgs inflation, the inllaton field is identified to the Higgs field h recently discovered at the 
Large Hadron Collider. But for being in agreement with data, the model requires that the Higgs is 
non-minimally coupled to gravity. This is the simplest model arising in the context of the standard 
model of particle physics. In the Einstein frame, the scalar field potential is of the form 

l/(0)=A4(^l-e-v^'^/“p>)^ (8.6) 

and the only parameter A is fixed by normalizing the scalar power spectrum amplitude to Planck 
measurement, so that there no remaining free parameter. The Higgs model is equivalent to the 
Starobinsky model, based on the action S = dx‘^ y/—g i^pi^ + R ^/). 


8.3 The best inflationary model after Planck+BICEP2 


A similar analysis has been performed for Planck plus BICEP2 [10], assuming that the B- 
modes detected by BICEP2 are of primordial origin, and using only the first four band powers (the 
others preferring a large value of the tensor-to-scalar ratio already strongly disfavored by Planck 
and WMAP). 

The principal finding is that there exist a net tension between Planck and BICEP2 data, most 
models favored by Planck alone being disfavored by BICEP2 alone (in terms of Bayesian evidence), 
and inversely. There nevertheless remains a series of models that are compatible both with Planck 
and BICEP2, and whose Bayes factor is reported in Pig. Note that the reference model is not 
Higgs inflation (which is not compatible with BICEP2) but a slow-roll model with £i, £2 and £3 as 
free parameters. In this case, it is found that large field models are favored and the simple potential 
V{(j)) = is one of the best models. 


9. Open problems 


In this section I briefly mention a discuss a few open issues related to inflation. 

9.1 Inflation vs. alternatives: 

Alternative to inflation exist and include string gas cosmology, matter bounces, ekpyrotic/cyclic 


scenarios (for a review, see e.g. Rev. [58]). But most of them do not solve simultaneously all the 
problems of the standard cosmological scenario, and thus the inflationary paradigm is often con¬ 
sidered as the best option, even if it must be seen as a model among others and in no way a theory 
already proved by observations. Even the detection of B-modes should not be considered as a 
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Displayed Evidences: 193 


Figure 8: Figure fom Q| displaying the Bayes factors for all inflationary scenarios listed in the ASPIC 
library []l5|], using Higgs inflation as a reference model. The bar color denote the class to which the model 
belongs (using Schwartz-Terrero-Escalante classification). Arrows indicate the maximum likelihood 


proof of inflation. Nevertheless, single field models of inflation have specific predictions such as 
the consistency relation between the scalar and tensor spectral indexes. Checking observationally 
that it is verified would strongly disfavor alternatives to inflation, and thus this will be one major 
objective of future CMB experiments. 

9.2 Reheating: 

At the end of inflation, the energy stored in scalar field potential must decay into standard 
model particles. It is thought that when inflation stops, the field start to oscillate around the min¬ 
imum of the potential. The coherent oscillations can be considered as a collection of independent 
scalar particles. If they couple to other particles, the inflaton can decay perturbatively to produce 
light particles. Another possibility is a phase of tachyonic preheating, which occurs in hybrid 
inflation scenarios due to the exponential grow of the modes of an auxiliary field during the wa¬ 
terfall phase. The modes become quickly non-linear and reheating occurs with the dissipation of 
interacting classical field oscillations 

As we have shown, the observable predictions for a single field model in slow-roll are easily 
calculated when one knows 4, the time at which the observable pivot scale exit the Hubble radius 
during inflation. But to find 4, one does not only needs to solve the dynamics of inflation but also 
all the subsequent expansion history. The simplest assumption is to consider that the reheating is 
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for Planck+BICEP2 and for models such that 7?.>1 
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Figure 9: Figure fom [ p^ displaying the Bayes factors for all inflationary scenarios listed in the ASPIC 
library [|^ compatible both with Planck and BICEP2, using a slow-roll model as a reference. As in the 
previous figure, the bar color denote the class to which the model belongs and arrows indicate the maximal 
likelihood value 


instantaneous and that a radiation era is immediately triggered at the end of inflation. But if this 
is not the case, then the reheating history affects the observable predictions for a given inflation 
model. One has therefore to derive reheating consistent constraints on its parameter space, e.g. by 
parametrizing the reheating phase by its duration and its mean equation of state (see e.g. ill])- 

9.3 Initial homogeneity problem: 

The question of how homogeneous must have been the initial conditions for inflation to be 
triggered has been tackled for more than twenty years by several authors 


63, 64,63]. 
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It has been shown for the simplest scalar field potentials that the Universe must have been initially 
homogeneous on scales larger than the Hubble radius. To draw this conclusion, relativistic simula¬ 
tions of the pre-inflation era have been conducted, either by using gradient expansion methods or 
by solving the full Einstein equations, in l-i-l dimensions (spherically symmetric case) or in 3-1-1 
dimensions. Thus inflation merely transforms one problem of homogeneity into another one, since 
an incredibly homogeneous initial state is required. This issue is called the Homogeneity Problem 
for inflation. 


9.4 Trans-Planckian problem: 


In inflation observable scales exit the Hubble radius about 60 e-folds before the end of infla¬ 
tion. But if inflation lasted for much more than 60 e-folds, then the observable wavelength modes 
were initially not only sub-Hubble but also sub-Planckian, in a regime where the physical princi¬ 
ples which underlie the calculations of the power spectrum are possibly not valid anymore. This 
occurs whereas the energy density itself is well below the Planck scale and so the homogeneous 
classical field dynamics is perfectly valid. There are several approaches to model trans-planckian 
effects, e.g. by modifying the dispersion relation, k —)■ w{k). One of them is to assume that Fourier 
modes are created when their wavelength equals some critical scale denoted by Me. This modifi¬ 
cation of the infrared behavior of the perturbation modes induce oscillations in the scalar power 
spectrum [F 


67] 


3i>^{k) = X I 


1 — 2|v|aocos 






(9.1) 


where Oq = H /M^, Ci, £2 the usual Hubble-flow parameters, visa complex number of modulus 
\x\ and phase (j) that parameterizes the initial conditions for the modes. Constraining the presence of 
oscillatory features in the scalar power spectrum therefore also constraints trans-Planckian physics. 


9.5 Avoiding the Big-Bang singularity 

Models of inflation alone are insufficient to describe the initial state of the Universe. They do 
not solve the initial singularity problem so that a quantum theory of gravitation is needed to describe 
the Universe at the Planck energy scale. But it is possible that the Universe was dominated by 
curvature just before inflation is triggered. In this case general relativity allows that the Universe 
has performed a classical bounce, thus avoiding the initial singularity, followed by a phase of 
inflation to explain the apparent flatness today. This scenario leads to specific signatures in the 
scalar power spectrum, taking the form of superimposed oscillations ||^]. Compared to the case 
of trans-Planckian effects, there are two important differences: for trans-planckian effects, the 
oscillatory term has a Ink dependance rather than a linear k dependance for bounding cosmologies. 

9.6 Eternal inflation: 

In most inflation models, there exist regions in the field space where the potential is so flat 
that its quantum fluctuations ~ H/In in a Hubble time dominate over its classical evolution 
^/H. Regions in the real space where quantum fluctuations push the field through values 
where the potential is more flat expand faster than others, and they can experience in turn quantum 
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fluctuations pushing the field towards more flat regions of the potential, which expand faster, and so 
on, so that they occupy a never-ending increasing part of the total volume. Inflation globally never 
ends, this is called the eternally self-reproducing regime. In this picture, our observable patch of the 
Universe emerges from an eternally inflating multiverse. The problem is that most models contains 
a self-reproducing regime. But eternal inflation leads to a collection of difficult problems like the 
measure non-normalisability and the apparent lost of unitarity (see e.g Refs. [1 ^ f/^ ] for a review 
on eternal inflation and related issues). 
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